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Abstract 

The two-loop 0(g 4 M 2 /M^) corrections are incorporated in the theoretical calculation of M, 
sm 2 9 w (M z ), and sin 2 6*^, as functions of M H . The analysis is carried out in a previously 
proposed MS formulation and two novel on-shell resummation schemes. It is found that the 
inclusion of the new effects sharply decreases the scheme and residual scale dependence of 
the calculations. QCD corrections are incorporated in two different approaches. Comparison 
with the world average of sin 2 6 l e e f pt f leads to M H = 127t^f GeV and M w = 80.367±0.048GeV, 
with small variations among the six calculations. 



* On leave from Dipartimento di Fisica, Universita di Padova, Padua, Italy. 



The theoretical calculation of the P^-boson mass M w , the MS parameter sin 2 6 w (M z ), 
and the effective mixing parameter sin 2 O^ff represents one of the most important applica- 
tions of the Standard Model (SM) at the level of its quantum corrections. We recall that 
sin 2 O^ff is determined from on-resonance observables at LEP and SLC, while sin 2 9 W (M Z ) 
is theoretically very important, particularly in the context of GUT studies. The aim of this 
paper is to present an accurate calculation of these parameters that, in the case of M w 
and sin 2 ^(M z ), includes the recently evaluated two-loop contributions of 0(g 4 M 2 /M^) [0]- 
The corresponding contributions for the evaluation of sin 2 9^fl are provided in the present 
paper. Although nominally subleading relative to the two-loop terms of 0{g A Mf /M^), these 
effects have been found to be of comparable magnitude [|]. Our strategy is to study the 
incorporation of the 0(g 4 M 2 /M^) contributions in different renormalization schemes, and 
to investigate the residual scale dependence. In fact, as it will be shown, the inclusion of 
the new results greatly decreases the scheme and residual scale dependence of important 
radiative corrections. QCD corrections are evaluated following two different approaches. 

We recall that, using the accurately known parameters a, G^, and M z as inputs, M w 
can be calculated from the expression |2j 

s 2 c 2 = A 2 /M 2 (1 - Ar), (1) 

where s 2 = 1 — c 2 is an abbreviation for sin 2 6 W = 1 — M^/Mj (M w and M z are the 
physical masses of the intermediate bosons), A 2 = na/y^G^ = (37.2802 GeV) 2 , and Ar is 
the relevant radiative correction ||. In order to incorporate the 0(g 4 M 2 /M^) effects, it is 
necessary to specify the renormalization scheme employed in the evaluation of the one-loop 
contributions. In Ref. M this was done in the MS framework of Ref. ||, where 

1 - Ar = [1 - A?v] s 2 /s 2 , (2) 
s 2 /s 2 = 1 + (c 2 /s 2 )(e 2 /s 2 ) Ap, (3) 
Ap = Re [A ww (Ml ) - c 2 A zz (M 2 )]_/M 2 w , (4) 

Ar w = -(25e/e) m +(e 2 /s 2 )f, (5) 
/ = [(ReA ww (Ml) - A ww (0))/Ml + V W + M 2 W B W ]_. (6) 

In these expressions MS denotes both the MS renormalization (i.e. the pole subtraction) and 
the choice [i = M z for the 't-Hooft scale, s 2 = 1 — c 2 is an abbreviation for sm 2 9 w (M z ), 2Se/e 
is the charge renormalization counterterm, e 2 = e 2 /(l + 2Se/e)jj^ the MS electromagnetic 
coupling at /i = M z , and A ww (q 2 ) and A zz (q 2 ) the transverse W and Z° self-energies ^-|J 
with g 2 = e 2 /s 2 factored out. In this paper we also factor out g 2 in the definition of Ap. As 
M t has been measured, it is simplest to use a pure MS subtraction, without decoupling the 
top contribution to 5e/e and the mixing self-energy A 1Z (M Z ) The term V w + M^B W in 
Eq. (H) represents vertex and box diagram contributions to ji decay, modulo a factor e 2 /s 2 [|]. 
QCD corrections to the self-energies are given, for example, in Ref. and are updated in 
this paper. 



1 



As shown in Ref. PJ, in the MS framework of Eqs. ([!]-[]) the 0(g 4 M 2 /M^) corrections 
are incorporated by adding two-loop irreducible contributions, Arff and {e 2 /s 2 )Ap^ 2 \ to 
the one-loop expressions Ar$ and (e 2 /s 2 )Ap^\ respectively. The correction (e 2 /s 2 )Ap^ 
also includes the previously evaluated irreducible contributions of 0(g 4 M^/M^) 0. 

In order to discuss the scheme dependence, we note that a resummation formula analogous 
to Eq. (0) can be obtained in the on-shell (OS) renormalization scheme of Ref. @. For brevity, 
we show how it follows from Eqs.([TJ-||). Combining Eqs. we have 



1 - Ar = 1 + (2Se/e) m + {c 2 / s 2 ) (e 2 /s 2 ) Ap - (e 2 /s 2 )f. 



(7) 



In order to express e 2 /s 2 in terms of G„, we note that Eqs. 



e 2 /s 2 = (G,/V2)8M ] 



(26e/e) 



MS 



) lead to 

(e 2 /s 2 )f 



Approximating e 2 / s 2 
Eq. (^), we have 



e 2 /s 2 in the square bracket and inserting Eq. 



1 - Ar 



1 + (2Se/e) m - -f 



1 + ^SMl^LAp 



V2 



(8) 

in the Ap term of 
(9) 



The leading M t contribution to Ap*- 1 ) is 3M 2 /(647r 2 M^), which is independent of s 2 . At 
the one-loop level (25e/e)y^ is also independent of s 2 and e [|, H, and (e 2 /s 2 )/W does 
not involve terms proportional to M 2 or large logarithms ln(M z /my), where rrif is a small 
fermion mass . However, non-leading one-loop contributions contained in (e 2 /s 2 )f^ and 
(G[j,/\/2) 8M^Ap^ do depend on s 2 . To obtain a resummation formula involving only on- 
shell parameters, we replace s 2 = s 2 [1 + (c 2 /s 2 ) (e 2 /s 2 ) Ap] in f(s 2 ) and Ap(s 2 ). Calling 
f{s 2 ) and Ap(s 2 ) the resulting functions, we obtain 



1 - Ar 



1 + 



25e 



MS 




- 3/> 



l + c 2 8M 2 .G, A _, 2; 



V2 



-Ap(s' 



(10) 



which is an on-shell counterpart of the MS expression of Eq. (§). We note from Eq. (||D that 
s 2 — s 2 = c 2 (e 2 /s 2 )Ap « 3c 2 x t + ... where x t = G ^M 2 / '8n 2 \J2 and the ellipses represent 
subleading contributions. Thus the replacement s 2 = s 2 [1 + (c 2 /s 2 ) 3x t + ...} in /^(s 2 ) and 
Ap^s 2 ) induces additional contributions of 0(g 4 M 2 /M|,). The corresponding functions, 
called f add and Ap a J d , respectively, are given in the Appendix. Therefore in the OS scheme 
we obtain 



where f^ 2 \s^ 



/ (2) (« 2 



Ap> 2 ) 

+ fadd( s2 ) 



/ 

Ap 



(2), 



/ 



+ Ap 



n( 2 )^ 2 ) 



(11) 

(12) 



Ap( 2 )(s 2 ) = Ap^(s 2 ) + Ap^is 2 ). The amplitude 



(2) 



(e 2 /s 2 )/( 2 ) is given by Eqs.(7a,b) of Ref. [| multiplied by (a/ns 2 )N c x t (N c = 3), while 
(G/t/y/2) 8M^ Ap^ 2 ^ is given by Eqs.(10a,b) of the same reference multiplied by N c x 2 . 
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An alternative OS resummation can be obtained by combining (e 2 /s 2 ) Ap in Eq. (|7|) with 
Eq. and once more replacing / — > f, Ap — > Ap. This leads to 

Ar = ArW + Ar^ + (^G M ^Ap(s 2 )J (1-Aa), (13) 

where Ar*- 1 - 1 is the original one-loop OS result of Refs. 0,f§, expressed in terms of a and 
s 2 , Ar^ = (e 2 /s 2 )f^ — (c 2 / 's 2 )(e 2 / 's 2 )Ap^ 2 \ and Aa is the renormalized photon vacuum 
polarization function at q 2 = M 2 . If in the last term we only retain the two-loop contributions 
proportional to M 2 and M 4 4 , Eq. ( |i~3"D reduces to 

Ar = Ar« + Ar» + ( £)' JV C *« - A^^T) . (14) 



The last terms in Eq. (13,14) represent higher order reducible contributions induced by re- 



summation of one-loop corrections, while Ar^ 2 -* contains the corresponding irreducible com- 
ponents. Eqs. (P| jroUIcSyi4T) satisfy the important property that, when inserted in Eq. ([[]), 
every factor of a in that equation is matched by a factor (1 + 25e/e)^ or, equivalently, 
(1-Aa)- 1 §. 

In the previous considerations we have set p = M z . However, one may also consider the 
case of a general p. Although the exact Ar, being a physical observable, is p- independent, 
we note that the resummation formulae of Eqs. (0,|TI],[nj) involve terms quadratic in gauge- 
independent, but p-dependent one-loop amplitudes. As a consequence, these resummation 
formulae contain a /i-dependence in 0(g 4 ). We have verified that the inclusion of the irre- 
ducible two-loop effects, and Ap^ in Eq. (0), and and Ap^ 2 ^ in Eqs. (|l(],[Uj), cancels 
the /^-dependence through 0(g 4 M 2 /M^), the order of validity of the calculation. There 
remains a significantly smaller /i-dependence of 0(g A ) without M 2 enhancement factors due 
to the fact that complete 0(g 4 ) corrections have not yet been evaluated. On the other hand, 
Eq. flUD is exactly /{-independent since it retains only the complete two-loop contributions 
proportional to M? and M 2 . We refer to Eqs.(|,|), Eq. (0), and Eq. flU) as the MS, OSI, 
and OSII schemes, respectively. Their numerical difference will give us a measure of the 
scheme dependence of Ar and the corresponding M w predictions. 

The effective parameter sin 2 9 l ^Jl is obtained from s 2 or s 2 by means of the relations 

sin 2 6 l e f f = k(M 2 ) sm 2 6 w (M z ) = k(M 2 ) s 2 (15) 

where k(q 2 ) and k{q 2 ) are the real parts of electroweak form factors evaluated at q 2 = M 2 , 
and s 2 and s 2 are related by Eq. (||). The amplitude k, evaluated in the MS scheme with 
p = M z , can be expressed as k = l + (e 2 /s 2 ) (Ak^ +Ak^). The one-loop contribution Ak^ 
is given in Ref . [[J , with the understanding that the top contribution to the mixing self-energy 
A 1Z (M 2 ) is not decoupled and, following the discussion of that work, certain two-loop effects 
induced by the imaginary parts of the self-energies are included. We recall that, for large M t , 
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AfcW grows like ln(M 2 /M|). The 0(g 4 M 2 /Ml) corrections are contained in (e 2 / s 2 ) Ak^\ 
which is given in Eq. ( |A3|) , multiplied by N c (a/4irs 2 ) 2 M 2 /(c 2 M 2 ). The on-shell amplitude 

(GJV2) 8 Ml \l - (e 2 /s 2 )/l (which follows 



k is obtained from k using Eq. (|3]), e 2 /s 2 
from Eq. (H)), approximating (e 2 / s 2 ) f (s 2 ) t 



8Mlf(s 2 ) in the square bracket, ex- 



pressing s 2 in AA;W(s 2 ) -f AA;^ 2 ^(s 2 ) in terms of s 2 and taking into account the additional 
0(g 4 M 2 /Ml) contribution Ak a 2 ] d induced by the latter shift. Neglecting two-loop terms 
without M 2 enhancement factors, this leads to 



1 + ^^Ak(s 2 ) 



1 | c 2 8M 2 w G, 
s 2 V2 



v/2 




Ap(s 2 



(16) 



where Ap is defined in Eq. (0) and AA;(s 2 ) = A£;« (s 2 ) + AA;( 2 ) (s 2 ) , AA;( 2 ) (s 2 ) = Ak^ (s 2 ) + 



Ak a 2 ] d (s 2 ). If we only retain two- loop effects with Mf and M 2 enhancement factors, Eq. (|T 



reduces to 
k 



1 + 



8MlG^ 



V2 



Ak(s 2 ) + - 2 Ap(s 2 ) + - 2 N c x t (AkV(s 2 ) - /«< 



;i7) 



which is exactly /^-independent. 

We briefly explain the strategies followed in the evaluation of the QCD corrections. When 
electroweak amplitudes proportional to the squared top-quark mass are expressed in terms of 
the pole mass M t , the coefficients of a n (t) (a = a s /n, n — 1,2) are quite large, afeature that is 
absent when they are expressed in terms of the running MS mass p, t = m t (p t ). Recent studies 
of the scale and scheme dependence of the leading M t -dependent electroweak corrections also 
suggest that the unknown QCD corrections to the 0(g A M 4 /Ml) contributions are significant 
when Mt is employed fllO| |. These observations suggest the strategy of parametrizing the 
electroweak amplitudes in terms of pt. We illustrate this approach in the most sensitive 
amplitude, namely the contributions to Ap and Ap proportional to powers of M t . In the 
case of Ap we have 



8M^ 
V2 



N, 



V2 8ti 2 



1+^ + ^^ 



8^2ti 2 



+ 



7= -tt, 

2V27T 2 



add 



+ 



where R add = 16n 2 Ap add / (N c x t ) and 



.MS 
°QCD 



-0.19325 a{n t ) -3.970 a 2 (p 



th 



(19) 



r H = M 2 /p 2 , r z = M 2 /p 2 , R(r H ,r z ) is the function given in Eqs.(10a,b) of Ref. U with 
M t — > pt and c 2 — > c 2 , and the ellipses stand for contributions not proportional to powers of 
p t . The coefficients in Eq. ( |T9"D are much smaller than in [[J 

Sqcd = -2.8599 a{M t ) - 14.594 a 2 (M t ), (20) 



the QCD correction when Ap is expressed in terms of M t . In the case of Ap we have an 
analogous formula to Eq. (|IB"D, except that one retains the s 2 dependence of R, R a dd is absent 
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and the couplings are expressed in terms of e 2 /s 2 . This strategy is extended to the other 
electroweak amplitudes using the perturbative relation M t = [l + (4/3) a(/i t ) + 8.236 a 2 (fj, t )] fh 
in the terms of 0(a) and 0(aa s ) of Ref. ||. As M t is important in the interpretation of the 
experiments, we employ 

Ht/M t = [1 + (8/3)a 2 (^) - 4.47a 3 (/i t )]/[l + (4/3)a(/i*) - 1.072a 2 (^*)], (21) 

where fl t = 0.252M t , //| = 0.0960M t . Eq. ([H]) avoids large coefficients and follows from 
the BLM optimizations of iXt/rht{M t ) and M t /rht(M t ) [fLl], with a(fl t ) and a(^) evaluated 
from a(M z ) using five active flavors. The corresponding FAC and PMS optimizations lead 
to nearly identical results [|ll[j. In this approach, fi t is obtained from M t via Eq. ( pT|) and 



inserted in the various amplitudes such as Eq. (TP8|). The second strategy followed in the 
paper is based on the conventional M t parametrization and Eq. (|20|) . In both approaches we 
include very small 0(aa 2 ) effects arising from the light isodoublets |l^] and from an inverse 



top mass expansion of the (t — b) isodoublet contribution [13| . As inputs in the Tables we 



employ M z = 91.1863 GeV, a s (M z ) = 0.118 0, Aa had = 0.0280 0, M t = 175 GeVft 



The 0(g 4 ) contribution to ReA zz (M 2 ) in Eq. (|) involving (lmA JZ (M 2 )) 2 and effects of light 
fermion masses in the one- loop corrections are also taken into account. 

Tables 1 and 2 compare the M w and sin 2 9 1 ^ predictions based on OSI (Eqs.( |Il||16| )), 
OSII (Eqs.flll 0)), MS (Eqs.(|,|,|l5D), with and without the complete 0(g 4 M 2 /Ml) cor- 
rections, respectively. Specifically, in Table 2 we neg lect the terms of 0(g 4 M 2 /Ml) in the 
functions / and Ap, as well as the contributions of the same order in the last terms of 
Eqs. ([14] ) and (|T7| ) and the second factor of Eq. (|I~6D. In both cases, the QCD corrections are 
evaluated using the /i r parametrization, explained in Eq. (^) et seq., and Eq. (|2~1"D. Table 
3 includes the 0(g 4 M 2 /M^) corrections but evaluates the QCD corrections using the M t - 
parametrization and Eq. (|20|) . In all cases the last column gives the values of sm 2 9 w (M z ) 
(MS evaluation). 

A number of striking theoretical features are apparent: (a) The incorporation of the 
irreducible 0(g A M 2 /M^) corrections greatly decreases the scheme dependence of the predic- 
tions, as measured by the differences encountered in OSI, OSII, and MS. For instance, the 
maximum differences in Table 2 are 11 MeV in M w and 2.1 x 10~ 4 in sin 2 9 l Sf, while in Table 
1 they are reduced to 2 MeV and 3 x 10~ 5 , respectively (similar very small differences can 
be observed in Table 3). (b) The differences between Table 2 and Table 1 are quite small in 
OSI. This means that this scheme absorbs a large fraction of the 0(g 4 M 2 /M^) corrections 
in the leading contributions, (c) Although the QCD approaches we have considered are quite 
different, Tables 1 and 3 show very close results. This is due to a curious cancellation of 
screening and anti-screening effects in the difference between the two formulations. 

Comparing the current world average sin 9 l e J} = 0.23165 ± 0.00024 with Tables 1 and 
3, and taking into account in quadrature the errors associated with Aa^ (±2.3 x 10~ 4 ) 

1 Results for different input values can be very easily obtained using a Mathematica package available 
from the authors. 
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and 5M t = ±6 GeV (±1.9 x 10 4 ), together with an estimate of missing higher order QCD 
corrections (±2.4 x 10~ 5 ) [I1J we obtain a determination of M H . As the six calculations 



of sin 2 in those tables are very close, it matters relatively little which one is employed. 
Averaging the six central values and choosing the errors to cover the range of the six cal- 
culations, we find M H = 127l^ 3 GeV, compatible with MSSM expectations. The strength 
of this determination rests on the fact that unaccounted two-loop effects are not enhanced 
by (M 2 /Ml) n (n = 1,2) factors or large logarithms ln(M z /mj).They are expected to affect 
sin 2 9^Fj in the fifth decimal, which is consistent with the variations observed in Tables 1 
and 3. Because of the approximate exponential dependence of M H on sin 2 this precision 
is important in order to obtain a reasonably accurate determination of M H . For instance, 
a 0.1% difference in the theoretical calculation of sin 2 9^1 induces a ~55% shift in the 
sin 2 9 1 ^ determination of M H and its la bound. Corresponding to the above determination 
of M H , one finds from Tables 1 and 3 the prediction M w = 80.367 ± 0.048GeV, in very good 
agreement with the current world average M w = 80.356 ± 0.125GeV. 
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hospitality during the summer of 1996, when part of this work was done, and M. Passera for 
useful calculations involving the resummation formulae. His work was supported in part by 
NSF grant PHY-93 13781. 



Appendix 

In the following we adopt the notation and conventions of Ref. |J. The additional contribu- 
tions to / and Ap, in units N c x t / (16n 2 ), are given by 

fl = y + 3^ + 4 c 2 B0[M 2 w , 0, M 2 W ] - + _L + 4 c 2 ) B0[M 2 w , M$, M 2 } 

(ll-8c 2 )lnc 2 /ll 1\ , A1 . 

+- tt^t hr + TT ln ^> ( A1 ) 



6 s 2 V 3 3 c 

Apfi = 5 -§ - ^ - ± I (l ± 26 c 2 ± 24 c 4 ) B0[M 2 z , M 2 W , M 2 W ] + 4 c 2 B0[M 2 w , 0, M 2 W , 

~ (y + 3^ + 4c 2 ) B0{M 2 w , Ml , M 2 ] - (1 + ^ - 8c 4 ) lnc 2 



, 2 38 34c 2 \ M? 2 3 - 62c 2 ±74c 4 ±36c 6 )\nzt 

+ cT 2 + In — f + — n " • A2 

1 3 3 J fi 2 9 c 2 



Fixing the one- loop result to the MS calculation of Ref. |§, the 0(g 4 M 2 /Ml) two-loop 
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contribution to AkS 2 \M z ) in units iV c /(i67r 2 ) (d/(47rs 2 )) M 2 /(c 2 M 2 ) reads [T§ 



-211 + 24 fa + 462 s 2 - 64 fa s 2 
432 



jj-ij 50[M|,M^,M^]-|lnc 2 
| (fa-4)y/fa(8g 2 -3)ff(fa) (6 + 27fa - 10 fa 2 + fa 3 ) (8 s 2 - 3) In fa 



108 



108 (fa - 4) 



1 2 A2 \ M 2 (3s 2 -2)ln^ , (fa-l)(8s 2 -3)0(f~ 
1 11 ~~ 2" + 18 - 



4 9 



+ 



18 (4 - fa) fa 



(A3) 



while the term to be added in the OS framework is given, in units iV c a;t/(167r 2 ), by 



Ak 



(2) 
add 



-238 c 2 



27 



+ 8c 4 - 2 c 2 V4c 2 - 1 (3 + 4 c 2 ) arctan( 



16 2 , 
— c lnzt 



+ 



4 c 2 



- 2 c 2 Wy(l) +4c>(c^) - 7c 2 In c l - — c 2 In 



Ml' 



(A4) 



where the functions fy{x) and ^y(x) are defined in Eqs.(6d) and (6e) of |17 
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M w (GeV) 


sin 2 effi 


"2 
S 




OSI 


OSII 


MS 


OSI 


OSII 


MS 


TV TO 

Mb 


65 


80.404 


80.404 


80.406 


0.23132 


0.23132 


0.23130 


0.23120 


100 


80.381 


80.381 


80.383 


0.23152 


0.23154 


0.23151 


0.23142 


300 


80.308 


80.307 


80.309 


0.23209 


0.23212 


0.23209 


0.23200 


600 


80.255 


80.254 


80.256 


0.23248 


0.23250 


0.23247 


0.23239 


1000 


80.216 


80.215 


80.216 


0.23275 


0.23277 


0.23275 


0.23267 



Table 1: Predicted values of M w , sin 2 0^/, and sin 2 6 w (M z ). OSI: Eqs.(0,0); OSII: 
Eqs.(|T3JT7D; MS: Eqs.(@,g,0). QCD corrections based on /i t -parametrization and Eq. Q2T]). 

M t = 175GeV, & S (M Z ) = 0.118, Aa had = 0.0280. 





M w (GeV) 


sin e eff 


s 2 


M H 


OSI 


OSII 


MS 


OSI 


OSII 


MS 


MS 


65 


80.409 


80.419 


80.417 


0.23132 


0.23111 


0.23124 


0.23114 


100 


80.385 


80.395 


80.393 


0.23154 


0.23135 


0.23145 


0.23135 


300 


80.311 


80.316 


80.318 


0.23212 


0.23203 


0.23203 


0.23195 


600 


80.257 


80.258 


80.263 


0.23251 


0.23247 


0.23242 


0.23235 


1000 


80.215 


80.214 


80.222 


0.23279 


0.23280 


0.23271 


0.23264 



Table 2: As in Table 1, but excluding 0(g 4 M 2 /M^) corrections (see text). 





M w (GeV) 


sm e eff 


s 2 


M H 


OSI 


OSII 


MS 


OSI 


OSII 


MS 


MS 


65 


80.405 


80.404 


80.406 


0.23132 


0.23134 


0.23130 


0.23121 


100 


80.382 


80.381 


80.383 


0.23153 


0.23155 


0.23152 


0.23142 


300 


80.308 


80.306 


80.308 


0.23210 


0.23214 


0.23210 


0.23201 


600 


80.254 


80.252 


80.254 


0.23249 


0.23252 


0.23249 


0.23241 


1000 


80.214 


80.213 


80.214 


0.23277 


0.23279 


0.23277 


0.23269 



Table 3: As in Table 1, but with QCD corrections based on M r parametrization and Eq. (^UP. 
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